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ABSTRACT 

This  report  presents  a  set  of  free-form  optimal  shapes  and  the  corresponding  stress 
concentration  factors  for  two  interacting  closely-spaced  holes  in  a  large  biaxially  loaded  plate. 
A  range  of  interaction  distances  and  two  distinct  biaxial  loading  cases  are  considered,  namely 
tensile  field  (remote  principal  stresses  have  the  same  sign)  or  mixed  field  (remote  principal 
stresses  have  the  opposite  sign).  The  optimal  shapes  are  obtained  using  the  finite  element 
analysis  based  gradientless  shape  optimisation  method  recently  developed  in  AVD.  In  a  key 
unexpected  result,  which  has  not  previously  been  identified,  it  has  been  discovered  that  the 
peak  stress  associated  with  the  optimal  holes  is  independent  of  the  interaction  distance,  and 
that  the  peak  stress  is  the  same  as  that  for  an  optimal  single  hole.  In  both  stress  field  cases,  the 
absolute  value  of  the  tangential  stress  is  piecewise  constant  around  the  hole  boundary. 
Compared  to  interacting  closely-spaced  circular  holes,  which  are  a  common  feature  in  aircraft 
structures  producing  significant  stress  concentrations  that  often  lead  to  premature  fatigue 
cracking,  the  optimal  hole  shapes  provide  for  very  large  reductions  in  peak  stress  (typically 
greater  than  30%),  which  results  in  a  substantial  increase  in  fatigue  life.  The  free-form  optimal 
shapes  are  presented  in  a  tabular  form  that  allows  them  to  be  used  readily  by  designers.  For 
the  tensile  biaxial  field,  the  optimal  shapes  are  smooth  and  non-circular.  For  the  mixed  biaxial 
field  case,  the  optimal  shapes  are  approximately  rectangular  with  sharp  comers  and  curved 
sides,  and  no  stress  singularities  are  present  at  the  comers. 
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Shape  Optimisation  of  Two  Closely-Spaced  Holes  for 

Fatigue  Life  Extension 

Executive  Summary 


In  recent  years.  Air  Vehicles  Division  (AVD)  has  been  undertaking  research  to 
improve  the  capability  to  optimise  the  precise  free-form  shapes  of  holes  and  fillets,  in 
the  context  of  life  extension  of  ageing  RAAF  aircraft  components.  Here  the  approach  is 
to  precisely  reshape  the  holes  to  an  optimal  shape  that  minimises  the  magnitude  of  the 
peak  stresses  on  the  hole  boundary,  and  thereby  inhibits  the  formation  of  fatigue 
cracks.  The  application  of  such  optimal  shapes  to  RAAF  aircraft  components  can 
provide  large  savings  in  the  cost  of  ownership,  through  an  extension  of  the  time 
intervals  between  costly  inspections.  For  example,  a  notable  recent  application  is  the 
shape  optimisation  of  holes  in  the  F-lll  wing  pivot  fitting.  While  there  has  been 
significant  progress  to  date  in  this  topic,  there  are  certain  technical  issues  that  need  to 
be  addressed  to  derive  the  most  benefit  from  the  shape  optimisation  approach.  One 
such  issue  is  the  effect  of  interacting  closely-spaced  holes,  which  are  a  very  common 
feature  in  RAAF  aircraft.  Typically,  standard  circular  holes  are  non-optimal  and  result 
in  significant  stress  concentrations,  which  can  lead  to  undesirable  fatigue  cracking  in 
components  early  in  their  service  life.  In  addition,  there  are  no  published  solutions  for 
optimal  hole  shapes  that  could  offer  a  practicable  alternative  to  the  use  of  circular 
holes. 

This  report  presents  practical  solutions  to  the  problem  of  determining  optimal  hole 
shapes  for  two  closely-spaced  holes  in  a  large  biaxially  loaded  plate.  During  this  work, 
enhancements  to  the  recently  developed  AVD  finite  element  based  gradientless  shape 
optimisation  procedure  have  been  implemented.  Two  distinct  biaxial  loading  cases  are 
considered,  namely  tensile  field  (remote  principal  stresses  have  the  same  sign)  or 
mixed  field  (remote  principal  stresses  have  the  opposite  sign).  The  resulting  free-form 
optimal  shapes  provide  very  large  reductions  in  stress  (typically  greater  than  30%),  as 
compared  to  the  initial  circular  holes.  For  the  tensile  field,  the  solution  shape  is  smooth 
and  non-circular,  and  the  stress  is  uniform  around  the  hole  boundary.  However,  for 
the  mixed  field  case,  the  final  optimised  shape  is  approximately  rectangular  with  sharp 
corners  and  curved  sides.  Here,  the  absolute  magnitude  of  the  tangential  stress  has 
been  rendered  piecewise  constant.  The  optimal  shapes  are  presented  in  a  tabular  form 
that  allows  them  to  be  used  readily  by  designers.  In  a  key  numerical  result,  which  has 
not  been  identified  previously,  it  was  found  that  the  optimal  shapes  eliminate  stress 
interaction  effects,  with  the  stress  concentration  produced  by  the  optimal  holes  for 
both  the  tensile  and  mixed  field  biaxial  loading  cases  being  identical  to  the 
corresponding  single  hole  optimal. 

The  procedures  developed  and  the  results  obtained  during  this  work  have  contributed 
to  an  enhanced  capability  in  AVD  for  determining  and  specifying  optimal  shapes  for 
extending  the  fatigue  life  of  reworked  holes  in  RAAF  aircraft  structural  components. 
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1.  Introduction 

Multiple  interacting  closely-spaced  holes  are  a  common  feature  in  aircraft  (e.g.  F-lll  wing 
pivot  fitting)  and  other  engineering  structures,  and  they  are  often  fatigue-prone  locations 
because  of  the  stress  concentrations  that  they  introduce.  As  a  result,  the  analysis  of  the  stress 
concentrations  caused  by  interacting  circular  holes  has  been  of  particular  interest  (Ling, 
1948;  Heywood,  1952,  1969;  Haddon,  1967;  Christiansen,  1968).  In  many  cases,  the  holes  are 
circular  in  shape  because  of  ease  of  manufacturing  considerations,  but  other  shapes  such  as 
ellipses  have  also  been  considered  to  try  to  reduce  stress  concentration  effects.  Ideally, 
designers  would  like  to  be  able  to  minimise  the  resultant  stress  concentrations  for  interacting 
holes,  and  the  associated  problem  of  determining  optimal  hole  shapes  has  important 
practical  ramifications  for  maximising  the  fatigue  life  of  aircraft  components.  As  a  result,  the 
problem  of  determining  optimal  shapes  for  interacting  holes  has  been  the  subject  of 
considerable  study. 

Cherepanov  (1974)  studied  two  interacting  holes  in  a  biaxially  loaded  plate  with  uniform 
loads  of  the  same  sign  and  a  constant  internal  uniform  pressure  p  acting  on  the  boundary  of 
the  holes.  The  concept  that  the  tangential  stress  was  constant  around  the  boundary  of  an 
optimal  hole  was  mentioned.  Although  an  analytical  solution  was  presented  for  what  was 
considered  at  the  time  to  be  the  case  of  closely-spaced  holes,  there  was  an  error  in  the 
expression  for  the  edge  separation  distance.  This  meant  that  in  practice  the  holes  were  not  at 
all  closely  spaced,  so  any  stress  interaction  effects  were  in  fact  minimal.  Vigdergauz  (1976) 
studied  plates  with  n  interacting  holes  in  a  biaxially  loaded  plate  with  uniform  loads,  with 
each  hole  loaded  by  a  uniform  pressure  p.  Although  example  solutions  of  optimal  hole 
contours  were  presented  for  n  =  2,  3,  4  and  6  holes,  the  loading  was  greatly  simplified  and 
consisted  only  of  the  internal  pressure  p  (no  remote  loading).  Because  the  presence  of  remote 
loads  is  such  a  typical  characteristic  of  aircraft  structures,  these  solutions  have  highly 
restricted  scope  for  any  possible  practical  applications.  A  significant  ancillary  outcome  of 
this  work  was  a  correction  to  the  expression  for  the  edge  separation  distance  in  the  hole 
geometry  equations  provided  previously  in  Cherepanov  (1974). 

Vigdergauz  (1977)  once  again  studied  a  biaxially  loaded  plate  with  uniform  loads  and  n 
holes,  with  each  hole  internally  loaded  by  a  uniform  pressure  p.  A  proof  was  given  that  the 
tangential  stress  around  an  optimal  hole  contour  is  a  minimum  compared  to  the  maximum 
obtained  on  any  other  hole  contours,  which  supported  the  assumption  that  had  been  utilised 
earlier  by  Cherapanov  (1974).  Example  optimal  hole  shapes  loaded  only  by  an  internal 
pressure  p  were  provided  for  n  =  3  and  4,  but  for  only  one  edge  separation  distance  per  case. 
Analytical  expressions  corresponding  to  the  optimal  hole  shapes  for  n  =  3  were  given. 
Banichuk  (1977)  also  studied  plates  with  n  holes,  providing  another  proof  that  an  equi- 
stressed  boundary  is  optimal,  but  no  hole  shapes  were  provided.  Vigdergauz  (1982) 
considered  the  case  of  two  interacting  holes  in  a  plate,  aligned  on  a  line  perpendicular  to  the 
direction  of  the  applied  load,  with  the  holes  themselves  also  loaded  by  a  constant  pressure  p. 
A  procedure  for  obtaining  an  approximate  numerical  solution  was  described  for  this  pair  of 
holes,  but  this  was  valid  only  for  quite  large  hole  separations,  where  any  interaction  effects 
would  be  small.  Some  shapes  and  stress  concentration  factors  as  a  function  of  the  interaction 
distance  were  provided  for  the  loading  case  where  the  internal  pressure  was  the  only 
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applied  load.  Vigdergauz  (1988)  again  revisited  the  problem  of  n  interacting  holes  in  a  plate, 
giving  optimality  conditions  but  no  optimal  shapes.  His  work  indicated  that  the  tangential 
stress  around  the  optimal  hole  boundary  is  the  same  for  all  cases,  i.e.  2  to  n  holes. 

Notwithstanding  all  of  the  above  activity  to  date,  no  paper  gives  the  stress  concentration 
factors  and  optimal  shapes  for  biaxially  loaded  plates  with  two  closely-spaced  interacting 
holes  for  a  range  of  interaction  distances.  To  this  end,  the  recently  developed  AVD 
gradientless  shape  optimisation  procedure,  which  is  based  on  (he  use  of  finite  element 
analysis  (FEA)  in  an  iterative  manner,  appears  to  be  ideally  suited  to  the  present  problem,  as 
it  has  previously  been  applied  successfully  to  problems  involving  fillets  and  single  holes. 

For  the  case  of  two  interacting  holes,  the  general  problem  geometry  and  the  corresponding 
notation  are  shown  in  Figure  1.  Here,  the  large  rectangular  plate  contains  two  centrally 
located  holes  of  height  h  and  width  zv.  For  circular  holes  of  radius  r  we  have  that  h  =  w  =  2r. 
The  plate  is  subjected  to  uniform  biaxial  remote  stresses  along  the  remote  edges  of  Si  in  the 
x-direction  and  S2  in  the  y-direction.  The  distance  between  the  edges  of  the  two  holes  is  e  and 
the  width  and  height  of  the  plate  are  24/z  and  2 7h,  respectively,  and  have  been  chosen  to  be 
much  greater  than  r. 

In  general,  the  local  stress  concentration  factor  Kt  around  the  boundary  of  a  stress 
concentrating  structural  feature  is  usually  defined  as 


K, 


*  local 


a, 


nominal 


(1) 


where  crlocal  is  the  value  of  the  local  stress  at  any  given  point  on  the  boundary  and  crnominal  is 
some  nominal  applied  stress.  It  is  possible  to  determine  a  peak  value  of  Kt  for  any  given 
configuration  of  stress  concentrator  geometry  and  loading  conditions. 

When  choosing  suitable  hole  configurations,  designers  often  rely  on  circular  or  elliptical  hole 
profiles,  for  which  the  maximum  stress  concentration  factors  are  provided  for  some 
standardised  geometries  and  loading  cases  (see  Pilkey,  1997).  However,  as  compared  to  the 
optimal  holes  obtained  in  the  present  work,  it  is  has  been  found  that  the  performance  of 
circular  or  elliptical  holes  is  poor  in  terms  of  the  associated  peak  value  of  Kt.  Nevertheless, 
(here  have  historically  been  a  number  of  reasons  why  circular  holes,  in  particular,  have 
continued  to  be  so  common:  (i)  ease  of  visualisation,  drawing  and  manufacture;  (ii)  there  are 
very  limited  analytical  solutions  for  interacting  true  optimal  shapes,  with  these  only 
available  for  restrictive  geometric  configurations  and/or  loading  cases  that  were  amenable 
to  analysis  (e.g.  see  Cherepanov,  1974,  and  Vigdergauz,  1982);  and  (iii)  accurate  numerical  or 
experimental  methods  are  needed  to  determine  optimal  hole  profiles  for  realistic  cases 
(however,  unfortunately,  no  such  transferable  data  exist).  In  relation  to  (iii),  it  appears  that 
only  recently  have  high-quality  finite  element  (FE)  based  shape  optimisation  procedures 
become  commonly  available.  The  fact  that  the  optimal  hole  profiles  need  to  be  reproduced 
accurately  in  order  to  deliver  the  full  optimal  Kt  value  is  not  considered  to  be  a  problem,  as 
modern  machining  techniques  can  readily  meet  the  precision  requirements. 

Optimal  hole  shapes  are  relevant  at  the  design  stage,  as  well  as  having  applications  in  local 
shape  reworking  during  service  to  remove  cracks  at  previously  non-optimal  holes.  For 
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example,  the  application  of  free-form  shape  optimisation  (for  holes  and  fillets)  in  a 
systematic  and  automated  manner  has  been  of  particular  recent  interest  in  AVD  in  order  to 
obtain  the  most  benefit  from  local  shape  reworking.  Examples  of  this  work  include  Kaye  and 
Heller  (1997),  Heller  et  al.  (1999),  Chen  et  al.  (2000),  Kaye  and  Heller  (2000),  Heller  et  al. 
(2000a,  2000b),  Waldman  et  al.  (2000),  Burchill  et  al.  (2002),  and  Heller  et  al.  (2002a,  2002b).  It 
is  important  to  note,  as  will  be  discussed  in  subsequent  sections,  that  an  essential 
requirement  for  optimality  is  achievement  of  uniform  tangential  stress  along  that  segment  of 
the  boundary  that  is  being  optimised  (see  Baud,  1934;  Heywood,  1945;  Neuber,  1972; 
Cherepanov,  1974;  Banichuk,  1977;  Schnack,  1979;  Vigdergauz  and  Cherkayev,  1986; 
Mattheck  and  Burkhardt,  1990;  Kaye  and  Heller,  1997;  Waldman  et  al.,  2001).  In  the  context 
of  reworking  existing  structural  components,  the  procedure  must  be  implemented  in 
accordance  with  the  constraint  that  reshaping  must  involve  material  removal.  However,  the 
procedure  could  equally  well  allow  for  material  addition,  which  would  be  appropriate  for 
shape  optimisation  at  the  design  stage. 

In  this  paper,  accurate  solutions  are  determined  for  the  optimal  hole  shapes  for  two  closely- 
spaced  holes  in  a  large  biaxially  loaded  plate.  Initially,  in  Section  2  the  key  features  of 
optimal  holes  and  notches  are  discussed,  along  with  a  description  of  a  recently  developed 
finite  element  based  gradientless  shape  optimisation  procedure,  which  is  employed  to  obtain 
the  optimal  shapes.  This  is  followed  in  Sections  3  and  4  with  detailed  stress  concentration 
results  for  holes  under  two  distinct  biaxial  loading  cases,  as  well  as  the  normalised 
coordinates  for  the  optimal  profiles.  These  two  biaxial  loading  cases  are  termed  tensile  field 
when  the  remote  principal  stresses  are  of  the  same  sign,  and  mixed  field  when  they  are  of 
opposite  sign.  A  range  of  interaction  distances,  e/h,  is  considered  for  both  the  tension  and 
mixed  field  cases. 


2.  Shape  optimisation  method  for  holes  and  notches 

2.1  General  features  of  optimal  holes  and  notches 

It  is  clear  that  a  key  feature  of  the  best  iterative  numerical  shape  optimisation  methods  is  the 
determination  of  a  boundary  shape  so  that  there  is  a  constant,  or  near  constant,  distribution 
of  tangential  stress  around  the  boundary.  The  existence  of  this  feature  has  been  shown 
experimentally  by  Baud  (1934)  and  Heywood  (1945),  and  numerically  by  Schnack  (1979), 
Mattheck  and  Burkhardt  (1990),  Heller  et  al.  (1999),  Waldman  et  al.  (2000),  Waldman  et  al. 
(2002)  and  Heller  et  al.  (2002a,  2002b),  to  be  consistent  with  the  aim  of  minimising  the  peak 
stresses.  Typically,  if  geometry  constraints  are  placed  on  the  possible  shape,  then  the  optimal 
shape  will  result  in  tangential  stresses  that  can  only  be  rendered  uniform  along  part  of  the 
boundary.  This  feature  appears  to  have  been  postulated  first  by  Neuber  (1972),  and  has  been 
used  successfully  in  shape  optimisation  of  fuel  flow  vent  holes  and  stiffener  runouts  in  the 
F-lll  wing  pivot  fitting  (Heller  et  al.,  2001a,  2001b)  to  increase  fatigue  life  by  reducing  the 
magnitude  of  stress  concentrations.  In  such  cases,  there  is  usually  a  sharp  comer  in  the 
optimal  geometry  at  the  end  of  the  uniform  stress  region  (but  without  stress  singularity), 
after  which  the  stress  decays  rapidly.  It  is  noted  that  those  analytical  solutions  that  are 
available  for  optimal  holes  typically  rely  on  special  geometric  or  loading  conditions  in  order 
to  be  able  to  obtain  solutions  where  the  tangential  stress  around  the  hole  is  rendered 
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constant  or  piecewise  constant  for  the  entire  hole  boundary  (e.g.  Cherepanov,  1974; 
Vigdergauz,  1982;  Vigdergauz  and  Cherkayev,  1986;  Cherkaev  et  al.,  1998).  For  this  reason, 
the  practical  applicability  of  such  analytical  solutions  is  judged  to  be  effectively  non-existent 
when  it  comes  to  their  potential  use  in  typical  aircraft  structures,  particularly  when  the 
applied  loading  is  complex. 

2.2  Shape  optimisation  using  gradientless  finite  element  method 

In  recent  work  in  AVD,  a  highly  effective  iterative  2D  gradientless  free-form  shape 
optimisation  algorithm  has  been  developed  and  implemented  using  the  standard  PAFEC 
Level  8.5  FEA  code  (see  Kaye  and  Heller,  1997;  Heller  et  al.,  1999;  Heller  et  al.,  2001a;  Heller 
et  al.,  2001b;  and  Waldman  et  al.,  2001).  The  PAFEC  FEA  code  is  invoked  directly  by  an  in- 
house  developed  Fortran  program  during  each  iteration,  and  it  has  been  used  for  all  the  FEA 
computations  described  in  this  paper.  The  general  approach  that  is  followed  here  is  based  on 
an  analogy  with  the  growth  behaviour  of  biological  structures,  where  material  is  added  at 
regions  of  high  stress  and/or  removed  where  the  stress  is  low,  a  concept  that  was  first 
identified  by  Thompson  (1917).  This  FE-based  approach  can  achieve  a  free-form  optimal 
shape  that  produces  a  constant  boundary  stress  distribution.  Apart  from  the  present 
approach,  shape  optimisation  (albeit  of  differing  effectiveness)  based  on  biological  growth 
analogies  has  been  utilised  experimentally  by  Heywood  (1945)  and  numerically  by  Mattheck 
and  Burkhardt  (1990).  However,  in  contrast  with  the  numerical  work  of  Mattheck  and 
Burkhardt  (1990),  the  present  approach  is  significantly  simpler,  both  in  terms  of  the 
approach  taken  for  defining  local  material  removal/  addition  on  the  free  boundary,  as  well 
as  in  terms  of  the  method  used  for  remeshing  between  iterations  in  order  to  avoid  mesh 
distortion.  Additionally,  a  particular  distinction  is  that  the  present  work  focuses  on  the  more 
difficult  aspect  of  constrained  reworking,  as  opposed  to  the  shape  optimisation  of  initial 
designs  that  may  be  less  constrained. 

In  the  context  of  hole  shape  optimisation,  the  present  gradientless  method  can  be  explained 
as  follows.  Consider  Figure  2a,  which  shows  a  plate  with  a  hole  (stress  concentrator)  with  a 
free  boundary,  T,  on  which  there  are  a  number  of  nodes  i  =  l,...,k,  under  an  arbitrary 
remote  loading.  For  an  arbitrarily  shaped  (non-optimised)  hole,  the  tangential  stress  will 
generally  vary  around  the  boundary,  as  shown  schematically  in  Figure  2b.  Here,  the 
tangential  stress  is  defined  to  be  the  non-zero  principal  stress  at  the  free  boundary.  For 
simplicity  here,  let  us  assume  that  all  the  boundary  stresses  are  positive  in  sign.  The  amount 
of  material  to  be  added  or  removed  at  any  point  on  the  boundary  is  taken  to  be  directly 
proportional  to  the  difference  between  the  local  tangential  stress  and  a  suitable  reference 
value  (Heller  et  al.,  1999).  This  process  is  then  repeated  iteratively  until  the  boundary  hoop 
stress  is  constant,  or  nearly  constant,  to  within  a  prescribed  tolerance.  The  amount  to  move 
any  given  node  i  on  the  boundary  at  each  iteration  (in  the  direction  of  the  outward  normal) 
is  specified  by  the  following  equation 

(2) 

^th 
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where  a  positive  value  of  di  indicates  material  addition,  crf  is  the  tangential  stress  at  node  i 
on  the  boundary,  crth  is  a  non-zero  threshold  boundary  hoop  stress,  c  is  an  arbitrary 
characteristic  length  such  as  the  initial  hole  radius,  and  s  is  an  arbitrary  step-size  scaling 
factor  typically  in  the  range  of  0.001-0.05  (see  Kaye  and  Heller,  1997;  Heller  et  al.,  1999, 
Heller  et  al.,  2001a;  Heller  et  al.,  2001b;  and  Waldman  et  al.,  2001).  The  parameters  c  and  s 
remain  fixed  throughout  the  shape  optimisation  process. 

In  relation  to  the  selection  of  the  parameter  ath ,  a  number  of  interesting  features  become 
apparent.  Here,  the  specific  magnitude  of  the  threshold  stress  will  define  the  type  of  shape 
optimisation  process:  (i)  if  we  select  <J th  to  be  less  than  the  maximum  stress,  but  more  than 
the  minimum  stress,  material  will  be  both  added  and  removed  at  various  locations  around 
the  boundary;  (ii)  if  we  select  crth  to  be  the  same  or  greater  than  the  maximum  stress,  only 
material  removal  will  occur;  and  (iii)  if  we  select  <7  th  to  be  the  same  or  less  than  the 
minimum  stress  only  material  addition  will  occur.  If  we  are  reworking  a  stress  concentrator 
in  an  existing  component,  we  must  use  option  (ii)  or  option  (i)  with  a  further  constraint  on 
boundary  movements  so  that  material  is  not  added  (i.e.  di  ^  0).  In  practical  problems  where 
there  are  geometrical  constraints  restricting  the  possible  shape,  the  optimal  shape  may  have 
two  or  more  constant  stress  regions.  For  some  problems  it  may  be  desirable  to  separately 
optimise  the  peak  stresses  in  two  different  regions,  and  this  can  be  accomplished  by 
choosing  two  different  values  of  a th  and  hence  optimising  each  region  independently,  albeit 
simultaneously.  This  approach  was  used  in  Heller  et  al.  (2001b).  It  is  also  noted  that,  by  the 
inclusion  of  absolute  value  signs  in  Eq.  (2),  the  general  algorithm  also  works  well  to 
minimise  compressive  peaks  (if  they  exist). 

The  shape  optimisation  method  described  above  is  highly  suitable  for  implementation  with 
any  standard  FEA  code  that  permits  the  movement  of  nodes  after  the  completion  of  an 
analysis.  For  such  an  iterative  procedure,  a  useful  parameter  that  can  be  used  to  monitor 
solution  convergence  is  the  normalised  range  in  Kt  around  the  hole  boundary.  This  can  be 
written  as 


where  Ktm“,  and  K*vg  are  the  maximum,  minimum,  and  average  values  of  Kt  around 
the  hole  boundary.  Clearly,  as  W  approaches  zero,  the  stress  becomes  more  uniform  around 
the  boundary  and  the  optimal  hole  shape  is  approached. 


2.3  Procedure  for  local  element  remeshing  and  hole  shape  specification 

In  the  present  work,  the  implementation  of  the  shape  optimisation  algorithm  with  the 
standard  FEA  code  PAFEC  was  used,  as  has  been  the  case  in  prior  work.  Using  a  valuable 
and  convenient  mesh  generation  feature  named  PAFBLOCKs,  only  the  nodal  movements  at 
title  stress  concentrator  boundary  need  to  be  specified.  For  two-dimensional  FE  analyses,  the 
mesh  generation  blocks  are  usually  four-sided  subregions  of  the  FE  model  consisting  of  a 
number  of  plate  elements.  Figure  3  shows  two  mesh  generation  blocks,  which  are  labelled  as 
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B1  and  B2.  Using  two  parameters,  Ni  and  N2,  each  mesh  generation  block  can  be  subdivided 
into  a  number  of  smaller  elements  in  the  spanwise  and  chordwise  directions  as  shown  in  the 
diagram.  The  parameters  Ni  and  N2  define  how  many  elements  are  to  be  created  within  each 
mesh  generation  block,  as  well  as  the  spacing  of  those  elements.  In  order  for  proportionately 
sized  internal  elements  to  be  created,  only  the  locations  of  the  corner  nodes  and  the  number 
of  elements  along  the  sides  of  the  blocks  need  to  be  specified.  The  moving  boundary  during 
the  shape  optimisation  is  denoted  by  the  grey  line,  while  the  edges  of  the  elements  within 
each  mesh  generation  block  are  denoted  by  the  dotted  lines.  For  each  mesh  generation  block, 
note  that  the  edge  of  the  block  that  lies  opposite  to  the  moving  boundary  remains  fixed,  and 
the  fixed  boundary  is  denoted  by  the  thick  solid  black  line  in  the  diagram. 

Within  the  PAFEC  code,  the  edges  of  each  mesh  generation  block  are  defined  by  two  corner 
nodes  and  a  midside  node  in  such  a  way  that  these  three  nodes  He  on  a  circular  arc  or  a 
straight  line  segment.  The  former  applies  for  the  initial  FE  mesh  consisting  of  a  circular  hole. 
These  multiple  circular  arcs  are  used  by  PAFEC  for  the  generation  of  intermediate  nodes.  At 
every  iteration  during  the  shape  optimisation,  the  nodes  on  the  free  boundary  that  define  the 
edges  of  each  of  the  mesh  generation  blocks  ,  are  completely  free  to  move  to  any  location  in 
accordance  with  Eq.  (2).  The  internal  elements  within  each  of  the  mesh  generation  blocks 
simply  shrink  or  grow  in  proportion  to  this  nodal  movement,  which  helps  to  maintain  mesh 
integrity  by  avoiding  mesh  distortion,  something  that  often  causes  problems  for  other  FE- 
based  shape  optimisation  programs.  To  some  extent,  the  selection  of  the  number  of  segments 
to  be  used  to  subdivide  and  model  the  moving  boundary  is  arbitrary.  However,  if  there  are 
not  enough  segments  then  this  will  be  readily  apparent  in  the  solution  stresses,  which  will 
not  be  uniform  along  the  hole  boundary  as  required  for  a  true  optimal  solution. 

Consider  the  case  where  there  are  n  element  comer  nodes  arbitrarily  spaced  along  a  typical 
side  of  a  mesh  generation  block  that  consists  of  a  straight  edge.  At  the  m*  iteration,  the 
coordinates  ( xjty  j )  of  the ;th  node  are  as  shown  in  Figure  3.  The  predicted  movement  of  the 

boundary  node  (x,,y,)  is  dx  -  (Ax,,  Ay,).  Then  the  new  coordinates  of  the  n  element  corner 
nodes  at  iteration  ra+1  can  be  computed  using 


*;+1  -  xj  +  a*; 


(  y."  _  y-m  \ 

A j 

xm  -x' 

yn  j 


y“+1=y;  +  Ayr 


yr-y : 


(4) 

(5) 


2.4  General  modelling  details 

For  the  analysis  of  the  example  problems  described  in  Section  3,  there  were  two  elements  per 
mesh  generation  block  along  the  moving  boundary.  However,  in  Section  4,  because  of  the 
presence  of  sharp  corners  in  the  optimal  shape  for  the  mixed  field  case,  only  one  element  per 
mesh  generation  block  was  used  along  the  moving  boundary. 
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All  of  the  FE  analyses  that  are  reported  on  here  were  undertaken  using  a  Hewlett-Packard 
K260  Series  9000  computer.  All  of  the  FE  meshes  consisted  of  8-noded  quadrilateral 
isoparametric  elements  and  linear  elastic  plane  stress  conditions  were  assumed,  with 
Young's  modulus  E  =  73  GPa  and  Poisson's  ratio  v  =  0.32  (which  are  material  constants  that 
are  typical  of  an  aluminium  alloy). 

2.5  Loadpath  visualisation 

An  important  issue  related  to  computational  efficiency  and  the  avoidance  of  mesh  distortion 
during  the  course  of  the  shape  optimisation  is  the  judicious  choice  of  the  initial  hole  shape. 
Here  it  is  desirable  to  have  an  indication  of  what  the  eventual  optimal  shape  will  be,  in  order 
to  guide  mesh  definition.  In  this  context,  based  on  some  prior  work,  loadpath  visualisation 
could  be  a  useful  tool.  A  novel  procedure  for  computing  and  visualising  2D  loadpath 
orientations  and  loadpaths,  using  FEA  results,  has  recently  been  developed  by  Kelly  and 
Elsley  (1993, 1995)  and  Kelly  and  Tosh  (1998).  Kelly  and  Elsley  defined  a  loadpath  to  be  the 
trajectory  taken  within  a  structure  by  a  unit  quantity  of  applied  load,  beginning  at  a  point  of 
application  and  ending  at  an  equilibrating  reaction.  They  introduced  the  concept  of  a 
hypothetical  force  "stream  tube",  which  is  depicted  in  Figure  4.  For  any  selected  resolved 
force  direction,  depicted  as  the  x-direction  in  Figure  4,  there  is  no  flow  of  load  across  the 
boundary  of  the  force  tube,  and  the  loadpath  is  therefore  bounded  by  lines  along  which 
there  is  no  contribution  to  the  force  in  the  x-direction.  Note  that  in  Figure  4  the  applied  and 
reaction  loads  Fxa  and  Fxb  axe  aligned  with  the  x-direction,  but  an  equivalent  diagram  can 
also  be  produced  for  the  forces  in  the  y-direction  (or  any  other  chosen  direction). 

In  more  recent  work  by  Waldman  et  al.  (1999,  2002b),  simple  explicit  equations  for  the 
computation  of  the  x-direction  and  y-direction  loadpath  orientations,  ax  and  ay,  at  any 
point  on  the  structure  have  been  derived,  and  these  are  as  follows: 

tana,=^  (6) 

ox 

tana  =—  00 

where  <JX  is  the  normal  stress  in  the  x-direction,  cry  is  the  normal  stress  in  the  y-direction, 

and  rxy  is  the  shearing  stress.  The  above  equations  have  been  applied  to  computing 

loadpaths  for  unconstrained  optimal  and  initial  non-optimal  single  holes  in  plates  under 
uniaxial  and  biaxial  loads  (Waldman  et  al.,  1999,  2002b).  Particular  features  of  these  optimal 
shapes  for  single  holes  were  identified,  namely  the  absence  of  loadpath  recirculation  regions 
in  the  plate.  The  loadpath  visualisation  methodology  developed  previously  is  used  here  to 
study  loadpaths  for  the  present  case  of  two  closely-spaced  interacting  holes.  It  is  considered 
desirable  to  do  so  in  order  to  assess  the  behaviour  of  loadpaths  when  multiple  holes  are 
involved,  as  this  particular  configuration  has  never  been  investigated  in  this  manner. 
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3.  Analysis  of  two  closely-spaced  holes  -  tensile  field 

3.1  Analysis  of  initial  geometry 

For  the  present  numerical  analyses,  the  remote  biaxial  stresses  were  both  tensile  with  Si  =  Si, 
and  the  hole  radius  was  taken  as  r  =  20  mm  ( h  =  40  mm).  Due  to  symmetry,  only  one  quarter 
of  the  plate  was  modelled  in  the  FEA.  The  FE  mesh  to  be  used  in  the  shape  optimisation 
consisted  of  2024  elements  and  it  is  shown  in  Figure  5a.  A  more  detailed  view  of  the  region 
around  one  of  the  holes  is  shown  in  Figure  5b,  including  the  node  numbering  along  the 
moving  boundary  and  the  distribution  of  mesh  generation  blocks  used  to  define  the  hole 
boundary.  The  mesh  generation  blocks  are  denoted  by  the  thicker  black  lines,  while  the 
mesh  generation  block  edges  that  are  fixed  are  denoted  by  the  thickest  black  line.  The  Kt 
value  around  the  hole  boundary  for  this  geometry  and  loading  condition  was  defined  as 

K,=j-  (8) 

where  cr.  is  the  value  of  the  local  tangential  stress  at  node  i  on  the  hole  boundary.  For  this 
case  with  an  edge  separation  distance  e/h  =  0.5,  the  peak  Kt  for  the  initial  geometry  was 
2.906  at  the  location  ( x ,  y)  =  (0,  -/i/4).  It  should  be  noted  that,  for  a  single  circular  hole  in  an 
infinite  plate  under  this  loading,  a  constant  value  of  Kt  =  2.0  occurs  around  the  entire  hole 
boundary. 

3.2  Shape  optimisation  analyses  and  results 

For  convenience,  the  shape  optimisation  was  undertaken  for  the  context  of  reworking, 
where  only  material  removal  is  allowed.  Thus,  at  each  iteration,  all  nodal  movements  on  the 
hole  boundary  were  constrained  such  that  di  <■  0.  The  constraint  on  these  nodal  movements 
was  achieved  by  selecting  the  threshold  stress  used  in  Eq.  (2)  to  be  equal  to  the  maximum 
tangential  stress  occurring  at  the  nodes  on  the  hole  boundary  at  each  iteration: 

<*&  =  Maxfo)  (9) 

1=1 

The  characteristic  length  c  was  defined  to  be  equal  to  the  radius  of  the  initial  circular  hole, 
while  the  step-size  scaling  factor  was  set  to  s  =  0.02.  The  other  constraint  on  the  shape 
optimisation  domain  required  that  the  minimum  allowable  distance  between  the  holes  be 
fixed,  which  for  the  V^symmetry  FE  model  means  that  y  ^  -e/2  for  any  point  on  the  hole 
boundary. 

A  typical  solution  geometry  is  shown  in  Figure  6,  and  it  was  obtained  after  374  iterations 
using  an  initial  geometry  consisting  of  two  circular  holes  with  edge  separation  distance 
e/h-  0.5.  Note  that  the  edge  separation  distance  between  the  two  optimal  holes  has  been 
reduced  to  e/h  =  0.4568  because  of  an  increase  in  the  height  of  the  holes  that  results  from 
material  removal  during  the  shape  optimisation  process  (the  width  of  the  holes  also 
increased,  from  an  initial  h/w  - 1.0  to  h/w  =  0.7182).  A  comparison  of  the  stress  concentration 
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factor  Kt  around  the  hole  boundary  for  the  optimal  shape  and  the  corresponding  circular 
hole  with  e/h  =  0.4568  is  presented  in  Figure  7.  For  the  optimal  shape,  the  Kt  values  at  the 
nodes  along  the  hole  boundary  have  been  reduced  to  an  essentially  constant  level,  and  they 
lie  within  a  range  of  2.013  to  2.009  (the  range  being  equivalent  to  0.36%  of  S2).  This 
represents  a  reduction  of  32.7%  compared  to  the  peak  stress  concentration  value  of  Kt  =  2.990 
for  the  two  interacting  circular  holes  with  e/h  =  0.4568  (determined  in  a  separate  analysis).  It 
is  clear  that  the  Kt  value  for  tire  two-hole  optimal  approaches  the  Kt  value  for  the  single-hole 
optimal  in  an  infinite  plate  for  this  biaxial  loading  case,  which  corresponds  to  a  circular  hole 
with  Kt  =  2.0. 

Figure  8a  shows  a  contour  plot  of  the  normalised  bulk  stress  distribution  around  the  initial 
circular  hole,  while  Figure  8b  shows  the  equivalent  results  for  the  optimal  hole.  Here,  the 
normalised  bulk  stress  K  is  defined  as 


A 

K  = 


K 

Sl+S2 


(10) 


Here  the  2D  bulk  stress  k  (which  is  the  same  as  the  first  stress  invariant)  is  defined  as 


K  —  CT x  +  (X y  —  (Tjj  +  <X22 


(11) 


where  crn  and  <x22  are  the  principal  stresses. 

A  comparison  of  Figures  8a  and  8b  reveals  a  large  difference  in  the  two  normalised  bulk 
stress  distributions,  with  the  bulk  stress  distribution  for  the  optimal  solution  demonstrating 
an  important  feature.  Within  a  tolerance  ranging  between  -0.5%  and  +1.4%  of  the  remote 
normalised  bulk  stress,  the  normalised  bulk  stress  contours  shown  in  Figure  8b  are 
essentially  uniform  for  the  entire  plate,  being  almost  identical  (i.e.  only  slightly  higher)  to 
those  for  the  plate  without  any  holes.  (A  separate  FEA  was  undertaken  for  the  no-hole  case.) 
This  is  the  requirement,  as  proposed  by  Bjorkman  and  Richards  (1976),  for  an  optimal  single 
hole  to  be  regarded  as  being  harmonic  for  loading  conditions  where  S1S2  >  0,  whereupon  the 
introduction  of  the  optimal  hole  must  leave  the  original  stress  field  in  the  absence  of  the  hole 
everywhere  unchanged.  The  fact  that  this  property  of  being  harmonic  also  applies  to  a  case 
with  two  interacting  holes  has  never  previously  been  explicitly  identified.  It  is  worth  noting 
that  the  optimal  shape  looks  similar  to  the  deloid  given  by  Bjorkman  and  Richards  (1979)  for 
the  single  hole  configuration;  however,  in  that  particular  case,  the  deloid  was  computed  for  a 
quite  different  load  case  involving  a  non-uniform  stress  field  with  linear  gradients.  It  is  clear 
that  the  abovementioned  requirement  for  optimal  harmonic  hole  shapes  will  only  be  exactly 
satisfied  for  infinite  plates.  This  is  because  finite-width  plates  will  always  have  a  higher 
average  bulk  stress  resulting  from  the  stiffness  loss  due  to  the  presence  of  the  hole. 

Using  the  present  method,  optimal  shapes  and  their  associated  Kt  values  were  then 
determined  for  three  additional  edge  separation  distances,  namely  e/h  =  0.2101,  0.9756  and 
1.4815.  For  these  cases,  the  edge  separation  distances  for  the  initial  circular  holes  were 
e/h  =  0.25, 1.0  and  1.5,  respectively.  The  Kt  results  for  the  complete  set  of  four  optimal  shapes 
are  given  in  Table  1,  together  with  the  Kt  values  that  were  calculated  for  the  equivalent  non- 
optimal  circular  holes.  The  results  for  two  circular  holes  for  a  number  of  additional  values  of 
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e/h  are  also  tabulated.  It  should  be  noted  that  the  peak  Kt  for  the  circular  holes  always 
occurred  at  the  location  corresponding  to  node  60  in  Figure  5b,  which  lies  on  the  symmetry 
line  x  =  0.  From  the  results,  it  is  evident  that  the  optimal  Kt  values  represent  a  large  reduction 
as  compared  to  the  corresponding  Kt  values  for  non-optimal  interacting  circular  holes,  with 
the  optimal  Kt  values  lying  in  the  range  2.015  to  2.010.  These  values  are  equal  to  the  average 
of  the  Kt  values  at  the  nodes  around  the  boundary  of  the  hole,  and  the  variation  was  less 
than  0.2%  of  each  average  value. 

Table  1  also  presents  a  comparison  of  the  Kt  values  that  were  determined  using  an  analytical 
solution  for  the  stress  in  an  infinite  plate  with  two  equal  circular  holes  (Ling,  1948)  and  the 
present  FEA  computations  for  various  e/h  values.  The  present  FEA  results  are  all  in  excellent 
agreement  with  those  determined  by  Ling  (1948),  although  in  each  case  the  Kt  value  from  the 
FEA  is  slightly  higher  than  the  analytical,  but  never  by  more  than  0.79%.  It  is  considered  that 
these  slight  differences  are  due  to  the  influence  of  the  finite  plate  width  on  the  FEA. 

Figure  9  shows  a  plot  of  Kt  versus  e/h  for  the  complete  set  of  results  for  the  optimal  holes 
and  the  non-optimal  circular  holes.  The  optimal  Kt  values  are  essentially  equal  to  2  (to  within 
0.75%),  with  the  slight  increase  above  the  value  of  2  being  due  to  the  plate  finite- width 
effects,  which  appear  to  increase  a  little  as  the  edge  separation  distance  is  decreased. 
Originally,  a  plate  one-third  the  size  of  the  one  shown  in  Figure  5  was  used,  but  this 
produced  a  Kt  value  that  was  approximately  5%  greater  than  2,  so  the  use  of  the  larger  plate 
was  settled  on. 

Hence,  the  key  result  obtained  here  is  that  optimal  shapes  that  eliminate  stress  interaction 
effects  between  two  closely-spaced  interacting  holes  have  been  determined  numerically. 
This  work  has  led  to  the  unexpected  finding  that  the  peak  Kt  for  two  interacting  optimal 
holes  is  identical  to  that  for  the  corresponding  optimal  shape  for  a  single  hole.  This  is  an 
important  feature,  which  does  not  appear  to  have  been  previously  noted  anywhere  in  the 
available  literature.  For  illustrative  purposes,  the  optimal  shapes  determined  in  the  present 
work  are  plotted  in  Figure  10,  where  it  can  be  seen  that,  as  the  edge  separation  distance 
decreases,  the  shape  of  the  two  separate  optimal  holes  approaches  that  of  a  single  circular 
hole  circumscribing  them  both.  For  the  benefit  of  transferability  to  designers,  die  normalised 
coordinates  for  the  family  of  optimal  shapes  are  given  in  Table  2.  Prior  work  by  Waldman  et 
al.  (2000)  has  shown  that  a  compact  function  cannot  represent  such  optimal  shapes 
accurately,  and  that  designers  should  instead  use  cubic  splines  to  interpolate  these 
coordinates. 

3.3  Comparison  with  analytical  solution 

It  is  noted  that  Cherepanov  (1974)  gives  optimal  solutions  for  plates  containing  two 
interacting  holes,  where  the  exterior  of  the  two  holes  in  the  Z  =  ( x ,  y)  plane  is  mapped 
conformally  onto  the  exterior  of  two  slits  lying  along  the  real  axis  of  the  complex  f  =  f  +  ir/ 
plane.  These  two  holes  are  aligned  horizontally  along  the  x-axis,  which  is  perpendicular  to 
the  layout  shown  in  Figure  1.  In  the  general  formulation  for  the  symmetric  case  under 
consideration  here,  the  width  of  each  slit  was  set  to  unity  and  the  two  slits  were  assumed  to 
He  within  the  bounds  1  <  f  <  2  and  -  2  <  f  <  -1 .  Under  the  conformal  mapping,  f  =  ±1 
corresponds  to  the  locations  of  the  hole  edges  that  are  closest  to  the  axis  of  symmetry,  while 
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^  —  +2  corresponds  to  the  locations  of  the  hole  edges  that  are  furthest  from  the  axis  of 
symmetry. 

However,  there  is  an  error  in  the  expression  for  Xb  in  Eq.  3.17  in  Cherepanov  (1974),  which 
was  discovered  and  subsequently  corrected  by  Vigdergauz  (1976).  The  complete  and 
corrected  equation  for  the  hole  shape  is  given  below,  where  c  —  2xb  is  the  interaction  (edge 
separation)  distance  between  the  two  holes,  p  is  a  constant  applied  normal  stress  around  the 
edge  of  each  hole,  and  a  homogeneous  field  of  constant  stress  acts  at  infinity,  where 

<7X  =  cr“  and  ay  =  cr"  . 


b^2E(<p,6Q°)  +  ^d2F(<p,6  0‘) -W(4-*2X#2 -1) 
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<p  =  arcsin 
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(12) 


Here  the  functions  F(<p,a)  and  E(<p,a)  are  the  elliptic  integrals  of  the  first  and  second 
kinds,  respectively.  The  constant  cx ,  which  gives  the  scale  in  tire  physical  x-y  plane,  is 
considered  to  be  real  and  positive,  and  remains  undetermined  by  the  formulation  of  the 
problem. 

For  the  particular  problem  being  considered  here,  we  have  p  =  0 .  The  edge  separation 
distance  is  simply  e  =  2x|#=1  =  2xB  and  the  height  of  the  hole  is  h  =  x\(_2  -  x|f=1  •  This  leads  to 
the  following  expression  for  the  edge  separation  distance  ej  h: 


(13) 
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where  dl  is  a  constant  and  has  been  calculated  to  be  equal  to  1.22840.  From  Eqs.  (12)  and 
(13),  it  is  evident  that  both  the  optimal  hole  shape  and  the  edge  separation  distance  are 
solely  functions  of  one  parameter,  b/a.  Consequently,  the  key  point  here  is  that  there  is  no 
explicit  way  to  determine  the  optimal  hole  shape  for  a  specified  edge  separation  parameter 
e/h  and  load  parameter  b/a  for  a  given  practical  problem,  a  factor  that  significantly  limits 
the  usefulness  of  the  solution. 

In  Cherepanov  (1974),  the  edge  separation  distance  that  was  depicted  for  some  of  the 
optimal  holes  was  small  relative  to  the  height  of  each  hole.  However,  after  applying  the 
correction  to  the  expression  for  xb,  it  turns  out  that  the  edge  separation  distance  between  the 
holes  is  actually  relatively  large.  This  can  readily  be  seen  in  Figure  11,  where  we  consider  a 
range  of  load  cases  S2/S1,  with  the  values  of  e/h  ranging  from  a  high  of  4.91  to  a  low  of  2.15. 
Thus,  based  on  the  behaviour  of  the  Kt  for  two  interacting  circular  holes  presented  in  Figure 
9,  any  interaction  effects  are  expected  to  be  small  for  the  biaxial  loading  cases  that  were 
considered  by  Cherepanov  (1974),  as  the  values  of  e/h  that  were  used  there  were  all  greater 
than  2.  We  note  here  that  the  optimal  shape  for  a  single  hole  in  a  tensile  biaxial  stress  field  is 
an  ellipse,  whose  aspect  ratio  h/w  is  equal  to  the  biaxial  load  ratio  S2/S1,  producing  a 
uniform  stress  distribution  around  the  edge  of  the  hole  (Durelli  and  Murray,  1943).  The  hole 
shapes  as  determined  by  Cherepanov  (1974)  are  in  fact  approximately  elliptical,  with  their 
aspect  ratios  h/w  almost  equal  to  the  biaxial  load  ratio  S2/ Si,  as  would  be  expected  for  pairs 
of  holes  for  which  the  edge  separation  distance  was  large  enough  for  any  interaction  effects 
to  be  negligible  (i.e.  as  for  a  single  optimal  hole).  Hence,  it  turns  out  that  Cherepanov's 
solution  is  for  holes  that  do  not  interact  to  any  appreciable  degree. 

3.4  Comparison  of  loadpath  distributions 

Resolved  x-  and  y-direction  loadpaths  for  the  initial  and  the  optimal  configurations  were 
determined,  and  are  given  in  Figure  12.  For  the  initial  case  of  the  circular  hole,  depicted  in 
Figures  12a  and  12b,  we  see  that  compaction  of  the  resolved  x-  and  y-direction  loadpaths 
occurs  at  the  two  peak  stress  locations,  respectively.  Although  it  was  initially  expected  to  be 
present,  recirculation  is  not  apparent  in  any  of  these  loadpath  plots.  This  is  believed  to  result 
from  the  fact  that  the  nominal  and  optimal  shapes  are  relatively  similar.  It  is  also  possible 
that  a  significantly  more  refined  mesh  would  be  needed  to  enable  finer  details  to  be 
observed.  In  prior  work,  Waldman  et  al.  (1999)  found  that  recirculation  was  a  notable  feature 
of  non-optimal  geometries.  However,  in  that  work,  there  was  a  significant  shape  difference 
between  the  initial  and  the  optimal  holes.  For  the  present  investigation,  the  plotted 
loadpaths  show  no  clear  distinguishing  features  between  the  initial  and  the  optimal  shapes. 


4.  Analysis  of  two  closely-spaced  holes  -  mixed  field 

4.1  Analysis  of  initial  geometry 

For  this  particular  analysis,  the  biaxial  remote  stresses  along  the  plate  edges  were  taken  as 
Si  =  -V2S2  (i.e.  a  mixed  field  with  remote  stresses  of  opposite  sign  and  different  magnitude), 
and  the  initial  interaction  distance  was  taken  as  e/h  =  0.5.  Once  again,  due  to  symmetry  only 
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one-quarter  of  the  plate  was  modelled  in  the  FEA,  and  2024  elements  were  used  to  define  the 
mesh.  The  FE  mesh  for  the  plate  with  the  initial  circular  hole  is  shown  in  Figure  13a,  where 
the  mesh  refinement  in  the  lower  right-hand  comer  was  required  because  of  the  high 
applied  shear  stress  in  that  region  (since  the  magnitude  of  Si  is  of  the  same  order  as  S2).  The 
Kt  value  at  the  nodes  around  the  hole  boundary  for  this  geometry  and  loading  condition  was 
defined  as 


where  S2  is  the  dominant  load  component,  and  cr,.  is  the  value  of  the  local  tangential  stress 
at  node  i  on  the  hole  boundary.  The  peak  Kt  value  for  the  initial  geometry  was  equal  to  3.132 
at  the  location  (x,  y)  =  (0.4985fc,  -0.7892ft)  on  the  hole  boundary.  This  compares  with  a  peak 
value  of  Kt  =  3.5  for  the  case  of  a  similarly  loaded  infinite  plate  containing  a  single  circular 
hole,  so  it  is  apparent  that  a  small  degree  of  stress  alleviation  has  resulted  from  having  two 
circular  holes  in  close  proximity  to  each  other  because  of  shielding  effects. 


4.2  Shape  optimisation  results  and  discussion 

The  shape  optimisation  parameters  were  taken  as  c  =  20  mm  and  s  =  0.02.  For  these 
iterations,  in  order  to  minimise  both  the  tensile  and  the  compressive  peak  (which  would  be 
desirable  in  the  case  of  fully  reversed  cyclic  plate  loading),  the  following  value  was  taken  for 
the  stress  threshold  in  Eq.  (2): 


(15) 


Here  the  stress  threshold  ath  is  equal  to  the  average  value  of  the  absolute  stress  around  the 
hole  boundary.  This  is  useful  for  computational  robustness  in  situations  where  overshooting 
of  the  true  optimal  solution  is  identified,  as  it  allows  material  to  be  added  as  well  as 
removed  during  the  course  of  the  shape  optimisation.  The  fact  that  the  final  shape  may 
encroach  on  the  initial  one  is  something  that  is  allowable  when  carrying  out  shape 
optimisation  for  an  initial  design.  However,  in  the  event  that  a  requirement  exists  for  the 
optimal  shape  to  lie  outside  some  specified  hole  shape  (such  as  when  reworking  only  was 
allowed),  then  it  would  be  necessary  to  start  from  a  larger  initial  hole  shape. 

The  final  optimised  finite  element  mesh  after  341  iterations  is  shown  in  Figure  13b,  and  an 
enlarged  view  of  the  final  solution  geometry  is  shown  in  Figure  13c.  It  can  be  seen  that  the 
optimal  shape  has  sharp  comers  and  is  approximately  rectangular  but  with  curved  sides, 
with  a  height-to-width  ratio  h/w  =  1.3288.  As  expected,  and  fully  in  accordance  with 
Eq.  (15),  material  addition  as  well  as  material  removal  has  taken  place  around  the  hole 
boundary.  This  is  evident  when  comparing  Figures  13a  and  13b,  as  most  of  the  near-vertical 
side  of  the  rectangle-like  optimal  shape  has  moved  in  leftwards  towards  the  y-axis  relative  to 
the  initial  circular  shape.  The  edge  separation  distance  has  also  been  reduced  from  an  initial 
e/h  =  0.5  to  e/ft  =  0.3224  during  the  shape  optimisation.  A  comparison  of  the  stress 
distribution  around  the  hole  boundary  for  two  interacting  circular  holes  with  e/ft  =  0.3224 


13 


DSTO-RR-0253 


and  the  optimal  hole  shapes  is  given  in  Figure  14.  It  can  be  seen  that  the  absolute  value  of 
the  stress  concentration  along  the  optimal  hole  boundary  has  been  reduced  to  what  is  an 
essentially  piecewise-constant  level  of  Kt  =  1.9156.  As  compared  to  the  peak  Kt  value  of 
3.1195  for  two  interacting  circular  holes  with  e/h  =  0.3224,  the  level  of  stress  concentration 
has  been  reduced  by  38.6%,  which  is  an  excellent  improvement. 

The  above  result  has  a  number  of  key  similarities  to  that  of  optimal  single  holes  in  an  infinite 
plate  under  mixed  loading  where  S1S2  <  0,  for  which  Vigdergauz  and  Cherkayev  (1986)  have 
determined  Kt  values  using  a  semi-analytical  method.  However,  they  did  not  present  any 
calculated  optimal  hole  shapes.  They  noted  that  the  hole  shapes  are  almost  rectangular, 
having  curved  sides  with  corners  of  a  definite  included  angle,  with  the  ratio  of  the  sides  and 
the  included  angles  being  dependent  on  the  applied  loading.  They  also  noted  that:  (i)  a 
required  condition  for  optimality  was  that  the  absolute  magnitude  of  the  tangential  stress 
around  the  optimal  hole  boundary  should  be  constant,  with  there  being  a  discontinuous  sign 
change  in  the  stress  at  the  comer  points  without  any  singularities;  and  (ii)  the  lack  of  any 
stress  singularities  is  possible  because  of  the  corresponding  discontinuous  change  in  the 
local  unit  tangent  vector  at  the  comers.  These  characteristic  geometric  features  of  the  optimal 
shapes  and  associated  stress  conditions  have  also  been  identified  in  the  present  work  for  two 
interacting  holes  in  a  plate  under  mixed  loading,  including  the  presence  of  sharp  comers 
and  the  absence  of  corner  singularities.  These  results  are  typical  of  optimal  shapes  obtained 
for  mixed  field  loadings  such  as  the  one  considered  here,  or  for  cases  where  geometric 
constraints  are  active  (e.g.  optimal  fillets  such  as  those  considered  by  Waldman  et  al.,  2001, 
and  constrained  holes  such  as  those  considered  by  Kaye  and  Heller,  1997). 

For  comparison  purposes,  in  Figure  15  the  predictions  of  peak  Kt  for  a  single  optimal  hole 
given  by  Vigdergauz  and  Cherkayev  (1986)  are  plotted  for  various  values  of  the  biaxial  load 
ratio,  -S1/S2.  Here  they  are  compared  to  some  other  available  solutions,  including  the 
present  two-hole  optimal  solution  of  Kt  =  ±1.92  obtained  for  -Si/  S2  =  0.5.  A  single-hole 
optimal  solution  of  Kt  =  ±2.84  obtained  for  -Si/S2  =  1.0  using  the  present  method  is  also 
shown.  As  can  be  seen,  this  result  is  very  close  to  the  approximate  numerical  solution  of 
Kt  =  ±2.73  obtained  by  Vigdergauz  and  Cherkayev  (1986).  A  key  point  to  note  here  is  that  the 
optimal  two-hole  result  is  very  similar  to  that  for  the  case  of  a  single  optimal  hole,  indicating 
that  the  interaction  effects  between  the  two  holes  are  negligible,  a  result  that  has  not 
previously  been  available.  The  kink  that  is  located  near  -S1/S2  =  0.3  on  the  curve 
corresponding  to  the  single-hole  optimal  shape  appears  to  indicate  that  there  is  a 
typographical  error  in  the  data  published  in  Vigdergauz  and  Cherkayev  (1986).  Judging  by 
the  overall  trend  exhibited  by  the  data,  it  would  appear  that  the  point  corresponding  to  the 
value  -Si/S2  =  0.311  and  Kt  =  1.62  is  the  one  that  is  in  error.  For  further  comparison,  two 
additional  published  results  are  given  for  the  biaxial  load  ratio  case  -Si/ S2  =  1  for  similar  but 
sub-optimal  holes:  Pilkey  (1997)  gives  Kt  -  3.2  for  a  single  square  hole  (with  rounded 
corners),  and  Durelli  and  Rajaiah  (1979)  give  K  =  3.0  for  a  single  double-barrel  quasi-square- 
shaped  hole  (also  with  rounded  corners).  As  expected,  both  of  these  values  are  somewhat 
higher  than  the  optimal  value  of  Kt  =  2.73  given  by  Vigdergauz  and  Cherkayev  (1986). 

For  the  mixed  field  case  with  two  holes  and  a  biaxial  load  ratio  of  -Si/ S2  =  0.5,  a  number  of 
different  configurations  were  optimised  to  produce  optimal  holes  for  a  variety  of  edge 
separation  distances,  e/h.  Figure  16  shows  a  plot  of  Kt  versus  e/h  for  the  optimal  holes  as 
well  as  the  non-optimal  circular  holes.  The  Kt  values  for  the  optimal  holes  are  reduced  by 
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more  than  35%  compared  to  the  equivalent  circular  holes,  which  is  a  very  significant 
reduction.  Furthermore,  the  optimal  holes  produce  an  essentially  constant  Kt  that  is 
independent  of  e/h  for  the  range  of  interaction  distances  that  were  analysed,  while  the 
results  for  the  circular  holes  show  a  small  increase  in  Kt  with  increasing  e/h.  For  the  benefit 
of  transferability  to  designers,  the  normalised  coordinates  for  the  family  of  optimal  shapes 
are  given  in  Table  3. 

Figure  17  shows  normalised  bulk  stress  distributions  (see  Eq.  (10))  for  the  initial  circular  hole 
as  well  as  the  final  optimal  hole.  Inspection  of  Figures  17a  and  17b,  corresponding  to  the 
initial  and  the  optimal  solutions,  respectively,  reveals  that  there  is  a  significant  difference  in 
the  bulk  stress  distributions.  Note  that  the  peak  values  of  normalised  bulk  stress  are 
significantly  reduced  for  the  optimal  case.  However,  unlike  the  case  involving  the  biaxial 
tensile  field  loading  for  which  SiS2  >  0,  the  bulk  stress  distribution  is  not  uniform  for  the 
optimal  shape  for  the  mixed  field  loading  for  which  S1S2  <  0,  and  hence  the  hole  is  not  a 
harmonic  hole.  A  hole  is  regarded  as  being  harmonic  if  the  bulk  stress  distribution  in  the 
plate  without  a  hole  present  is  unaffected  by  the  inclusion  of  the  hole  (Bjorkman  and 
Richards,  1976).  Leading  on  from  this,  it  was  initially  envisaged  that,  for  the  mixed  field 
loading  for  which  SiS2<0,  the  difference  in  the  direct  stresses  (ax  -ay)  throughout  the 

plate  might  be  unaffected  by  the  presence  of  the  optimal  hole,  leading  to  a  new  class  of 
harmonic  holes.  However,  it  was  found  that  this  was  not  the  case.  Hence,  it  appears  that  it  is 
doubtful  that  any  further  work  is  required  to  extend  the  concept  of  optimal  harmonic  holes 
to  the  mixed  loading  case.  However,  the  present  work  confirms  numerically  that,  as  far  as 
optimality  is  concerned,  the  essential  requirement  is  piecewise  uniformity  of  the  tangential 
stress  along  the  edge  of  the  hole  boundary. 

4.3  Comparison  of  loadpath  distributions 

Resolved  x-  and  y-direction  loadpaths  for  the  initial  and  the  optimal  configurations  were 
determined,  and  are  shown  in  Figure  18.  For  the  case  of  the  initial  circular  hole.  Figure  18a 
reveals  that  compaction  of  the  resolved  r-direction  loadpaths  occurs  at  the  two  locations  of 
peak  compressive  stress  at  the  top  and  bottom  of  the  hole,  and  a  strong  recirculation  region 
is  evident  in  the  peak  tension  stress  region  on  the  right-hand  side  of  the  hole.  The  reverse  of 
this  behaviour  is  evident  for  the  resolved  y-direction  loadpaths  in  Figure  18b.  Conversely, 
for  die  optimal  shape.  Figure  18c  shows  no  recirculation  in  the  r-direction  loadpaths,  while 
recirculation  is  apparent  above  and  below  the  hole  for  the  y-direction  loadpaths,  and 
compaction  of  the  loadpaths  is  quite  evident  in  both  cases.  Comparing  Figures  18a  and  18b 
to  Figures  18c  and  18d,  the  recirculation  regions  are  significantly  less  pronounced  for  the 
optimal  shape  as  compared  to  the  initial  circular  shape.  This  is  consistent  with  the  findings 
by  Waldman  et  al.  (1999)  in  prior  work  for  single-hole  geometries.  It  is  believed  that  the 
recirculation  that  is  evident  for  this  optimal  shape  results  from  the  fact  that  sharp  comers  are 
present,  which  produce  a  rapid  change  in  sign  of  the  local  stress  field  that  may  act  like  a 
boundary  for  the  loadpaths. 
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5.  Conclusion 

In  this  report,  accurate  profiles  of  free-form  optimal  shapes  and  their  associated  stress 
concentration  factors  have  been  presented  for  large  plates  containing  two  closely-spaced 
interacting  holes.  Two  distinct  biaxial  loading  cases  and  a  range  of  hole  interaction  distances 
have  been  investigated.  In  a  key  unexpected  result,  which  has  not  previously  been 
identified,  it  has  been  discovered  that  the  peak  stress  associated  with  the  optimal  holes  is 
independent  of  the  interaction  distance,  and  that  the  peak  stress  is  the  same  as  that  for  an 
optimal  single  hole.  Overall,  it  has  been  demonstrated  numerically  that  the  nature  of  the 
distribution  of  tangential  stress  around  the  optimal  hole  boundaries  depends  critically  on 
the  sign  of  the  product  of  the  remote  biaxial  loads  Si  and  S2.  For  the  tensile  biaxial  loading 
cases  where  S1S2  >  0,  the  tangential  hoop  stress  is  constant  around  die  hole  boundary.  For 
the  mixed  field  biaxial  loading  cases  where  S1S2  <  0,  the  absolute  value  of  the  tangential 
hoop  stress  is  piecewise  constant  around  the  hole  boundary.  This  characteristic  behaviour 
for  die  mixed  biaxial  loading  case  was  first  identified  by  Vigdergauz  and  Cherkayev  (1986). 
Compared  to  interacting  closely-spaced  circular  holes,  which  are  a  common  feature  in 
aircraft  structures  producing  significant  stress  concentrations  that  often  lead  to  premature 
fatigue  cracking,  the  optimal  hole  shapes  provide  for  very  large  reductions  in  peak  stress 
(typically  greater  than  30%),  which  results  in  a  substantial  increase  in  fatigue  life. 

The  free-form  optimal  shapes  are  presented  in  a  tabular  form,  which  allows  them  to  be  used 
readily  by  designers.  For  the  tensile  biaxial  loading  case,  the  optimal  shape  is  non-circular, 
and  meets  the  definition  for  a  harmonic  hole  (i.e.  the  bulk  stress  field  in  the  loaded  plate 
without  a  hole  is  unchanged  by  the  inclusion  of  the  optimal  hole),  which  has  never  been 
identified  in  any  prior  work.  For  the  mixed  field  biaxial  loading  case,  the  final  optimised 
shape  is  approximately  rectangular  but  with  curved  sides,  and  with  sharp  comers.  When 
coupled  with  the  discontinuous  change  in  the  local  unit  tangent  vector  at  the  comers,  the 
discontinuous  nature  of  the  tangential  stress  either  side  of  the  corners  does  not  produce  any 
corner  stress  singularity  (Vigdergauz  and  Cherkayev,  1986). 

Loadpath  visualisations  for  the  tensile  biaxial  field  case,  where  S1S2  >  0,  found  that  no 
loadpath  recirculation  regions  were  present  for  either  the  circular  or  the  optimal  hole  cases. 
This  possibly  occurred  because  the  circular  and  optimal  shapes  could  be  considered  as  being 
relatively  similar  to  each  other.  For  the  mixed  field  biaxial  loading  case,  where  S1S2  <  0,  it 
was  found  that  recirculation  was  significantly  less  pronounced  for  the  optimal  shapes  as 
compared  to  the  initial  shape,  and  this  is  consistent  with  prior  work  for  unconstrained 
single-hole  geometries.  It  is  believed  that  the  existence  of  a  recirculation  region  for  the 
optimal  hole  shapes  for  S1S2  <  0  results  from  the  presence  of  sharp  comers,  which  produce  a 
rapid  change  of  sign  in  the  local  stress  field  that  may  act  like  a  boundary  for  the  loadpaths. 

The  procedures  that  were  developed  and  the  results  obtained  during  this  work  have 
contributed  to  an  enhanced  capability  in  AVD  for  specifying  optimal  shapes  for  extending 
the  fatigue  life  of  reworked  holes  in  RAAF  aircraft  structural  components. 
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Table  1:  Kt  values  for  two  interacting  holes  of  circular  and  optimal  shape  in  a  large  plate  under  remote 
biaxial  loading  Si  =  S2  for  various  e/h  values. 


e/h 

Kt  for  circular  holes 

Kt  for  optimal 
holes  from 
present  FEA 

Ling  (1948) 

Present  FEA 

Difference 
FEA  vs  Ling 

0.2101 

— 

3.955 

— 

2.015 

0.2500 

— 

3.682 

— 

— 

0.3500 

— 

3.259 

— 

— 

0.4568 

— 

2.990 

— 

2.011 

0.5000 

2.887 

2.906 

+0.65% 

— 

0.6000 

— 

2.755 

— 

— 

0.7500 

— 

2.595 

— 

— 

0.9756 

— 

2.437 

— 

2.010 

1.0000 

2.411 

2.426 

+0.62% 

— 

1.2500 

— 

2.322 

— 

— 

1.4815 

— 

2.258 

— 

2.010 

1.5000 

— 

2.254 

— 

— 

2.0000 

2.155 

2.172 

+0.79% 

— 

4.0000 

2.049 

— 

— 

—  - 

7.0000 

2.018 

— 

— 

— 

00 

2.000 

— 

— 

— 
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Table  2:  Coordinates  of  optimal  shapes  for  two  interacting  holes  in  a  large  plate 
under  remote  biaxial  loading  Si  =  S2  for  various  e/h  values. 


Si 

=  s2 

e/h  = 

0.2101 

e/h  = 

0.4568 

e/h  = 

0.9756 

e/h  = 

1.4815 

Kt  = 

2.0150 

Kt  = 

2.0110 

Kt  =  : 

2.0100 

Kt  =  : 

2.0100 

h/iv  - 

0.6015 

h/w  = 

0.7182 

h/w  = 

0.8601 

h/w  = 

0.9139 

x/h 

y/h 

x/h 

y/h 

x/h 

y/h 

x/h 

y/h 

0.0000 

0.1050 

0.0000 

0.2284 

0.0000 

0.4878 

0.0000 

0.7407 

0.0555 

0.1050 

0.0484 

0.2286 

0.0428 

0.4888 

0.0412 

0.7420 

0.1110 

0.1050 

0.0969 

0.2292 

0.0854 

0.4918 

0.0822 

0.7458 

0.1665 

0.1050 

0.1453 

0.2304 

0.1279 

0.4968 

0.1230 

0.7521 

0.2220 

0.1050 

0.1937 

0.2322 

0.1700 

0.5040 

0.1632 

0.7610 

0.2775 

0.1050 

0.2421 

0.2349 

0.2117 

0.5135 

0.2028 

0.7726 

0.3330 

0.1050 

0.2904 

0.2387 

0.2528 

0.5254 

0.2414 

0.7868 

0.3885 

0.1051 

0.3385 

0.2442 

0.2930 

0.5401 

0.2790 

0.8038 

0.4440 

0.1054 

0.3864 

0.2518 

0.3320 

0.5575 

0.3152 

0.8235 

0.4995 

0.1059 

0.4337 

0.2624 

0.3696 

0.5780 

0.3498 

0.8459 

0.5550 

0.1073 

0.4799 

0.2769 

0.4053 

0.6015 

0.3824 

0.8711 

0.6104 

0.1104 

0.5242 

0.2964 

0.4386 

0.6283 

0.4128 

0.8989 

0.6654 

0.1180 

0.5656 

0.3216 

0.4693 

0.6580 

0.4407 

0.9293 

0.7185 

0.1342 

0.6025 

0.3530 

0.4968 

0.6908 

0.4658 

0.9620 

0.7649 

0.1645 

0.6337 

0.3900 

0.5209 

0.7262 

0.4878 

0.9968 

0.7981 

0.2090 

0.6586 

0.4316 

0.5411 

0.7638 

0.5066 

1.0335 

0.8184 

0.2607 

0.6770 

0.4764 

0.5573 

0.8034 

0.5219 

1.0718 

0.8288 

0.3152 

0.6891 

0.5233 

0.5694 

0.8444 

0.5337 

1.1113 

0.8313 

0.3706 

0.6953 

0.5714 

0.5774 

0.8864 

0.5418 

1.1517 

0.8276 

0.4260 

0.6962 

0.6198 

0.5813 

0.9290 

0.5463 

1.1926 

0.8184 

0.4808 

0.6922 

0.6681 

0.5812 

0.9718 

0.5471 

1.2339 

0.8046 

0.5345 

0.6837 

0.7158 

0.5771 

1.0143 

0.5444 

1.2750 

0.7863 

0.5869 

0.6709 

0.7625 

0.5693 

1.0564 

0.5381 

1.3157 

0.7641 

0.6378 

0.6543 

0.8080 

0.5579 

1.0976 

0.5283 

1.3558 

0.7382 

0.6869 

0.6340 

0.8520 

0.5431 

1.1377 

0.5153 

1.3949 

0.7089 

0.7340 

0.6104 

0.8943 

0.5250 

1.1765 

0.4991 

1.4328 

0.6763 

0.7789 

0.5837 

0.9347 

0.5039 

1.2137 

0.4799 

1.4692 

0.6407 

0.8215 

0.5541 

0.9731 

0.4800 

1.2491 

0.4579 

1.5041 

0.6024 

0.8616 

0.5218 

1.0092 

0.4534 

1.2826 

0.4331 

1.5370 

0.5614 

0.8991 

0.4870 

1.0429 

0.4243 

1.3140 

0.4060 

1.5680 

0.5181 

0.9337 

0.4500 

1.0741 

0.3930 

1.3431 

0.3765 

1.5968 

0.4725 

0.9655 

0.4109 

1.1027 

0.3596 

1.3698 

0.3449 

1.6233 

0.4251 

0.9942 

0.3699 

1.1286 

0.3244 

1.3940 

0.3114 

1.6474 

0.3758 

1.0198 

0.3273 

1.1516 

0.2875 

1.4156 

0.2762 

1.6688 

0.3250 

1.0422 

0.2832 

1.1718 

0.2491 

1.4346 

0.2396 

1.6877 

0.2729 

1.0612 

0.2379 

1.1890 

0.2095 

1.4507 

0.2016 

1.7037 

0.2197 

1.0769 

0.1916 

1.2031 

0.1689 

1.4640 

0.1626 

1.7170 

0.1655 

1.0892 

0.1444 

1.2141 

0.1274 

1.4744 

0.1227 

1.7273 

0.1107 

1.0980 

0.0966 

1.2220 

0.0853 

1.4818 

0.0822 

1.7348 

0.0555 

1.1033 

0.0484 

1.2268 

0.0427 

1.4863 

0.0412 

1.7392 

0.0000 

1.1050 

0.0000 

1.2284 

0.0000 

1.4878 

0.0000 

1.7407 
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Table  3:  Coordinates  of  optimal  shapes  for  two  interacting  holes  under  remote  biaxial  loading 
Si  =  -V2S2  for  various  e/h  values. 


Si  =  -V2S2 


e/h  =  l 

0.2754 

e/h  =  0.3224 

e/h  =  0.3834 

e/h  =  0.6162 

e/h  =  0.9499 

e/h  =  1.3406 

Kt  =  1.9410 

Kt  =  1.9156 

Kt  =  1.9015 

Kt  =  1.9022 

Kt  =  1.9151 

Kt  =  1.9477 

hfw  = 

1.2235 

h/iv  = 

1.3288 

h/w  = 

1.5159 

h/w  = 

1.5497 

h/w  = 

1.7862 

h/w  = 

1.8398 

x/h 

y/h 

x/h 

y/h 

x/h 

y/h 

x/h 

y/h 

x/h 

y/h 

x/h 

y/h 

0.0000 

0.1377 

0.0000 

0.1612 

0.0000 

0.1916 

0.0000 

0.3081 

0.0000 

0.4750 

0.0000 

0.6703 

0.0413 

0.1382 

0.0424 

0.1621 

0.0416 

0.1925 

0.0406 

0.3088 

0.0396 

0.4755 

0.0374 

0.6706 

0.0825 

0.1397 

0.0848 

0.1647 

0.0831 

0.1953 

0.0811 

0.3110 

0.0792 

0.4769 

0.0748 

0.6715 

0.1237 

0.1425 

0.1270 

0.1689 

0.1244 

0.1997 

0.1215 

0.3144 
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Tangential  stress  at  boundary  node 
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1  i  k 

Nodal  position  around  contour 

(b) 


Fig.  2.  Explanation  of  moving  boundary  method  for  case  of  an  open  hole  in  a  remotely  loaded  plate, 
(a)  Geometry  and  notation  with  hole  boundary  shape  defined  by  k  nodes  on  contour  T .  (b)  Typical 
tangential  stress  distribution  around  hole  boundary  contour. 
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Fig.  3.  Use  of  automatic  mesh  generation  blocks  (B1  and  B2)  to  define  FE  mesh  for  shape  optimisation 
ofFEA  model. 
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Fig.  5.  Typical  finite  element  mesh  for  Vi-symmetry  model  of  large  rectangular  plate  containing  two 
closely-spaced  circular  holes,  (a)  Complete  finite  element  mesh,  (b)  Detail  of  mesh  around  circular 
hole  showing  mesh  generation  blocks  and  node  numbering  around  the  hole  boundary. 
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Fig.  6.  Detail  of  mesh  of  Vi-symmetry  FE  model  of  optimal  hole  shape  for  two  interacting  holes  under 
remote  tensile  biaxial  loading  Si  =  S2  with  final  e/h  =  0.4568. 
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Stress  cone 


Fig.  7.  Stress  concentration  factor  around  hole  boundary  of  Vi-symmetry  model  of  two  interacting 
holes  in  a  plate,  for  the  case  of  a  circular  hole  and  an  optimal  hole  for  e/h  —  0.4568  and  remote  tensile 
biaxial  loading  Si  =  S2. 
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(b) 


Fig.  8.  Contour  plots  of  normalised  bulk  stress  k  for  1 'A-symmetry  model  of  two  interacting  holes  in  a 
plate  for  e/h  =  0.4568  and  remote  tensile  biaxial  loading  Si  =  S2.  (a)  Non-optimal  circular  hole, 
(b)  Optimal  hole. 
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Minimum  edge  separation  distance,  e/h 


Fig.  9.  Variation  in  stress  concentration  factor  for  two  interacting  circular  and  optimal  holes  in  a 
large  plate  as  a  function  of  minimum  edge  separation  distance  between  the  two  holes,  e/h,  for  tensile 
biaxial  remote  loading  Si  =  S2. 


Fig.  12.  Loadpathsfor  \ A-symmetry  model  of  two  interacting  non-optimal  circular  holes  and  optimal 
hole  shapes  in  a  large  plate  for  remote  tensile  field  biaxial  loading  Si  =  S2  and  e/h  =  0.4568.  (a)  x- 
direction  and  (b)  y-direction  loadpathsfor  circular  holes,  (c)  x-direction  and  (d)  y-direction  loadpaths 
for  optimal  holes. 
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Fig.  13.  Finite  element  mesh  for  Vi-symmetry  model  of  Urge  rectangular  plate  containing  two  closely- 
spaced  holes  for  mixed  biaxial  field  remote  loading  Si  =  -V2S2.  (a)  Initial  circular  hole  with  e/h  =  0.5. 
(b)  Optimal  hole  with  e/h  =  0.3224.  (c)  Close-up  view  of  mesh  for  optimal  hole. 
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Stress  concentration  factor,  K 
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Biaxial  load  ratio,  -S-1/S2 


Fig.  15.  Comparison  of  stress  concentration  factors  for  single  optimal  holes  for  different  ratios  of 
biaxial  loading  (mixed  field),  including  results  for  two  optimal  interacting  holes  for  -S1/S2  =  1  and 
-$l/S2  =  V2. 
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Fig.  17.  Contours  of  normalised  bulk  stress  k  for  'A-symmetry  model  of  two  interacting  holes  in  a 
large  plate  for  mixed  biaxial  field  remote  loading  Si  =  -'ASi  and  e/h  =  0.3224.  (a)  Non-optimal 
circular  hole,  (b)  Optimal  hole. 
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Fig.  18.  Loadpaths  for  Vi-symmetry  models  of  two  interacting  non-optimal  circular  holes  and  optimal 
holes  in  a  large  plate  for  mixed  biaxial  field  remote  loading  Si  =  -V2S2  and  e/h  =  0.3224.  (a)  x- 
direction  loadpaths  and  (b)  y-direction  loadpaths  for  the  circular  hole,  (c)  x-direction  loadpaths  and 
(d)  y-direction  loadpaths  for  the  optimal  hole. 
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